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We present recent results in the charm sector from the CLEO Collaboration. The data used 
were collected by the CLEO II and III detectors at the Cornell Electron Storage Ring (CESR). 
Measurements of form factors and branching fractions in semileptonic decays of D-meson and 
charmed-baryons, searches for CP violation and rare decays, and Dalitz analyses of D-meson 
decays are discussed, as well as future plans for the CLEO-c experiment at CESR-c. 



1 Charm physics at CLEO 

Since the 1980's the primary goal of CLEO Collaboration was the study of T-resonances and 
£?-meson physics. The charm sector received steady attention as well. Among ~ 350 CLEO pub- 
lications about 18% were devoted to the decays of charmed mesons and about 14% to charmed 
baryons. At this time we continue to study all aspects of charm including semileptonic decays of 
mesons and baryons, Dalitz plot analysis of three-body D-meson decays, fragmentation, inclu- 
sive decay rates, doubly and singly Cabibbo suppressed decays, searches for rare and forbidden 
decays, CP-violation, etc. 

The data using in our analyses were collected with three detector configurations^^, 

CLEO LI Oct.1989 - Apr.1995 4.7 fb~ l , 
CLEO LL.V Nov.l99h - Fefr.1999 9 fb~ l , 
CLEO III JuZ.2000 - Mar.2003 ~ 16 fb' 1 , 

at Cornell e + e~ Storage Ring (CESR) with y/s ~ 10 GeV. At this time about half of the 
CLEO III data is available for analysis. 



2 Recent results overview 

2.1 First Observation of the Exclusive Decays — > A7r + 7r + 7r~7r° and A+ — ► Aw7r + 

This is a first publication^! based on CLEO III statistics, 4.1 fb _1 at T(4S*). The goal is to 
search for the inclusive multi-body baryonic decay A+ — > A7r + 7r + 7r _ 7r°. In spectator quark model 
this decay is equivalent to the D-meson decays observed with high rate, D + — > i^ 7r + 7r + 7r~7r° 
(BR = 5.4±?:2%), D° -> K-iT + iT + iT-iT {BR = 4.4 ± 0.4%). In this analysis we used the 
power of CLEO III reconstruction capabilities including K/p/ir particle identification using 
RICH &: dE/dx. We apply high momentum cut, P(A+) > 3.5 GeV/c, to suppress combinatoric 



background. Other cuts were quite standard. A clear signal from the decay A+ — > A7r + 7r + 7r~7r° 
is seen in the invariant A7r + 7r + 7r~7r mass spectrum, Fig. ^ 

A search for substructure in ir + ir~ir° combinations shows a signal from A+ —* Acj7t + , uj — > 
7r + 7T _ 7r° and a contribution from A+ — > A?77r + , r] — > 7r + 7r~7r° at ~ 3cr level. The later decay 
has previously been observed in the r\ — > 77 decay mode ^. We find that the decay A+ — > 
A7T + 7r + 7r _ 7r° is saturated by the kum + and A?77r + intermediate states. We find NO evidence for 
<jjix higher mass resonances. This study is limited by the allowed phase space. 
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Figure 1: Invariant A7r + 7r + 7r~7r° mass spectrum. 
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Figure 2: Definition of kinematic variables 
for the decay — > Ae + v e . 



We use the decay A^~ — > pK n + for systematic cross-checks and normalization, that gives 
the basic results for ratios of branchings 
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2.2 Form Factors and Search for CP Violation in the Decay A+ — ► Ke + u e 

The Heavy Quark Effective Theory (HQET) gives a reliable prediction on heavy to light quark 
transition for A- type baryons using the V + A hadronic current, 

< A 2 |< +A |A 1 >= u 2 [/ 1 ( (7 2 ) 7m (1 - 75 ) + f 2 (q 2 ) A 7m (1 - 75 )]m, 

parameterized in the frame of Korner &: Kramer (KK) model^. Thus the 4-fold differential 
rate of the decay A+ — > Ae + v e , 
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can be expressed as a function of four kinematic variables, q 2 , cos#a, cos 6\y , X ( see Fig. [^1 for 
explanation), and two parameters, pole mass, M po i e , and formfactors ratio, R = f2(l 2 )/ fiil 2 ), 
where fi{q 2 ) = { J£ q fy^ )2 (1 - q 2 max / 'M 2 ole ) 2 is taken in dipole form. 

To estimate the A+ momentum we use three sources of information: (i) direction of the 
event trust axis, (ii) the momentum balance equation P 2 + ~ (Pa + P e + P u ) 2 , (iii) weights from 
known fragmentation of A+. 

In order to extract R and M po i e we use 4-D maximum likelihood fit^and get R = -0.31 ± 
0.05 stat ± 0.04 syst and M pole = 2.13 ± 0.07 stat ± 0.10 syst GeV/c 2 . Note, that M pole from the 
fit is consistent with M D *+ value used as a parameter of KK model. The KK model predicts 
R = -0.25 at M po i e = M D ,+ =2.11 GeV/c 2 . 



We also measure the decay asymmetry parameter, a A + , which appears in parameterization of 

the differential rate after integration over x an d cos Ow that gives dq i^ os g A oc (l+a A +«A cos 6\), 
where a\ is a similar decay asymmetry parameter for A. We obtain from the fit a A + = — 0.85=L 
0.03 stat ± 0.02 syst @ < q 2 >= 0.67 {GeV/c) 2 for world average a A = 0.642 ± 0.013." 

The decay asymmetry can be obtained for two charge conjugate states separately, a A +a\ = 
—0.561 ± 0.026 s t a t, a A ~+a A = —0.535 ± 0.024 stat . These two values allow us to estimate a CP 

a A + JrCC + 

violation in this decay, A A + = c _ a _ c = 0. 01 =L0.03 s t a t ±0.013,^=1=0.02^. We find no evidence 

At A+ 

of CP violation in terms of the decay asymmetry parameter. 

2.3 First Search for Flavor Changing Neutral Current Decay D° — > 77 

Standard Model (SM) predicts the rate for the D° — » 77 decay of ~ 10~ 8 or less. Non SM 
extension, i.e. gluino exchange in SUSY, might enhance this rate by two orders of magnitude. 
A measurement of this rate is a good test of new physics beyond the SM. 

For this analysis^! we use statistics of CLEO II &: II. V. Combinatoric background is sup- 
pressed by the tagging process D* + — > D°ir + using the energy release variable Q = M(D* + ) — 
M(D°) — m n +, Fig. 12b. The decay D° — > 7r°7r° with measured rate is used for systematic 
cross-checks and normalization, Fig. The ratio of efficiency is determined from Monte Carlo. 
e(77)/e(7r 7r°) = 1.59 ± 0.05. We found 19.2 ± 9.3 signal events from the fit, and set an upper 
limit 

B(D° -» 7 7 )/B(D° -» tt°7t°) < 0.033, B(D° -» 77) < 2.9 x 10~ 5 @ 90%C.L. 
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Figure 3: Energy release in the decay 
D*+ D°it+ for (a) D° -> 7r°7r°; (b) 
D° 77. 



Figure 4: Dalitz plots for the decay D° -> Ksn + n~ for (a) D° , 
M RS = M{K° s -k-); (b) D°, Mrs = M(K° s n + ). 



24 Dalitz Plot Analyses of the Decays D -> PPP at CLEO 

The Dalitz plot analyses of D-meson decays in to three pseudoscalars, D — > PPP, are interesting 
for many reasons. 

• The Dalitz plot clearly shows the intermediate states and their interference. 

• These sorts of decays accounts for a large fraction of the total D meson decay rate. 

• It is a good place to study scalars, /o(980), ao(980). 

• In these decays there are indications on the not well established light scalar mesons, cr(500), 
k(400 - 600) (see E791, FOCUS, CLEO recent results). 



• An opportunity to search for CP violation in amplitudes and phases, that may be more 
visible than rate asymmetries. 

For the Dalitz plot analyses discussed here we use the CLEO II. V 9 fb _1 data sample, the large 
sample which has well understood resolutions and systematics. We also use a "generic" (all 
known decays included) and signal Monte-Carlo simulation in order to estimate an efficiency 
variation across the Dalitz plot. 

We exploit D* — > Dit s i ow tagging process in order to suppress background and to define an 
initial flavor of D° or £)°-meson. We use other standard for CLEO event selection algorithms, 
and ~ 2 — 3a cuts on energy release variable, Q = M(D*) — M(D) — M n , and invariant D-meson 
mass, M(D). 

The Dalitz plot is a scatter plot of two (of three) invariant masses squared, i.e. versus 
mf 2 - A general formalism of the matrix element parameterization and unbinned maximum 
likelihood fit is published in detail in Ref. Thus for the matrix element we use a sum of 
amplitudes of non-resonant and all possible resonant terms, 

M = a nr e i(pnr + a r e ltfir Aj(123\r), 

resonances 

where a r e lipr are the complex parameters (fit parameters) and ^4j(123|r) is a spin dependent 
Breit -Wigner amplitude for appropriate intermediate resonance, r. 

Based on parameters extracted from the fit one may define for each resonance a Fit Fraction = 

^ P ' sl ^|A4Rp!p dPP ' and an Inte 9 rated CP Asymmetry, A C p = § PS ^]^^T\M~^ d ' DP / $PS dVP 
The modified complex parameters can be also used in order to search for CP violation in terms 
of amplitudes and phases, a r e %Vr — ► a r e %Sr ^1 ± ^e^j, where a r and 5 r - CP conserving 

parameters; b r and 4> r - CP violating parameters, showing the difference in amplitude and phase 
between D° and D°. 



2.5 Search for CP Violation in the Decay D° — > Kgir^n 

The results of the Dalitz plot analysis of D° — > KgiT + iT~ were published!*^ for clearly observed 
and measured 10 modes, K*~ (892)vr+ , ^(1430)-vr+, K 2 *(1430)-7r+, K*(1680)-vr+, K%p, K° s uj, 
X£/o(980), Kg/ 2 (1270), Kg/ Q (1370), and the "wrong sign" K* + (892)^-. 

In present consideration we use the same clean sample of 5299±73 signal events with ~ 2% 
of background to search for CP violation (see Fig. The SM predicts for this decay a CP 
violation in rate at the level of ~ 10~ 6 due to the K° — K° mixing. Hence, an observation of 
CP violation in D° — > KgTr + ir~ would be a strong evidence for non-SM processes. 

The results will be published soon in terms of parameters, a r , 5 r , b r /a r , (f> r , Fit Fraction and 
Interference Contribution for 10 intermediate states and non-resonant amplitude. Preliminarily, 
we find NO indication to CP asymmetry in integrated parameter, Acp = —0.39 ± 0.34, and in 
all other types of CP asymmetries. 



2.6 Dalitz Plot Analysis of the Decay D° — > tt + tt tt° 

The decay D° — > ir + ir~ir° is interesting for a couple of reasons: 

• SM predicts a large CP violation in rate ~ 10~ 3 that is close to our sensitivity; 

• E791 found in D + — > 7r + 7r + 7r~ decay strong evidence for cr(500) — > 7r + 7r~ in the inter- 
mediate state. We would expect to observe this state or its charged isospin partner in 
D° -> vr+vr-vr . 
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Figure 5: The Dalitz plot for D° -» Figure 6: The m*_ x+ Figure 7: The m^ +ir „ Figure 8: The m^^ 
7r + 7r~7r° decay. projection. projection. projection. 



Table 1: Results of the Fit for the decay D° — > tt + tv tt°. Statistical and systematic errors are shown respectively. 



Interm. state 


Amplitude 


Phase (°) 


Fit Fraction (%) 


p°TT° 

Non resonant 


1 (fixed) 
0.56±0.02±0.07 
0.65±0.03±0.04 
1.03±0.17±0.31 


(fixed) 
10±3±3 
-4±3±4 
77±8±11 


76.5±1.8±4.8 
23.9±1.8±4.6 
32.3±2.1±2.2 
2.7±0.9±1.7 



Distributions of Q and M(D) shows ~ 1100 signal events with ~ 18 ± 3% of background. 
The Dalitz plot and three projections are shown in Figs[5HEl The observed K^n contribution 
is consistent with previous measurements and is not allowed to interferee with the other con- 
tributions to the Dalitz plot. Preliminary results of the fit parameters are listed in Tabled 
Our conclusions: Acp = O.OI^q'q? (stat) db 0.09(syst) and there is NO CP asymmetry in Acp 
at the level of our sensitivity; there is NO statistically significant difference in amplitude and 
phase between D°/D° — ► 7r + 7r _ 7r°; non-resonant contribution is small; there is NO indication 
on cr(500) fit fraction at the level of ~1%. 



2.7 First Observation and Dalitz Plot Analysis of the Decay D° — ► Kg7]ir° 

In the decay of D° — ► Kgnir one would expect a strong manifestation of D° — ► Kgao(980) , if 
the scalar resonance ao(980) indeed contributes to the complementary D° — > Kg(K + K~) decay 
along with /o(980). There is no result in PDG for D° — > K^nn yet. 

To study this decay we use the mode: Kg — > tt + tt~, rj — ► 77, it — ► 77. The decay 
D° — > Kgit is used for systematic cross-checks and normalization. We observe a clean signal 
in Q (Fig. EJ) and M(K s t]tt°) (Fig. Il()|) distributions and have preliminary measured the ratio of 
branchings 

BR(D° -> KSr/Tr ) 



In spite of large background ~ 30% we are going to perform a Dalitz plot analysis for the 
decay D° — > K^nn (see Fig. II 1|) . Our preliminary observations confirms that the iCgao(980) is 
a dominant intermediate state, the K* (892)77 * s seen with a significant interference term. To get 
a good fit to the observed data another contribution, such as non-resonant or high mass states, 
is needed. We are working on completing this analysis. 




3 CLEO-c and CESR-c project news 

The NSF has approved a five year program of charm physics called CLEO-c and CESR- 

(ID] The 

first year of this program has been successfully completed. Since October 2002 to March 2003 the 
crucial part of this project a single super-conducting (SC) wiggler had been successfully tested 
at CESR and five other wigglers have been assembled and prepared for installation. During a 
current shutdown since March 3 to July 2, 2003 the 6 SC wigglers will be installed in CESR. 
CLEO III will be upgraded to CLEO-c by replacing the silicon vertex detector (SVX) with a 
6-layers drift chamber with large stereo angle (ZD). 

The program of CLEO-c scans can be briefly listed as follow: 

• 2003 Act I: ^(3770) — 3 fb -1 ; 30M events, 6M tagged D (310 times MARK III). 

• 2004 Act II: yfi ~ 4100 MeV — 3 fb" 1 ; 1.5M D S D S , 0.3M tagged D s (480 x MARK III). 

• 2005 Act III: i/>(3100) — 1 fb -1 ; 1 Billion J/ip (170 times MARK III, 20 times BES II). 

This statistics is required for precise measurement of branching ratios, decay constants and other 
SM parameters in charm sector. 

4 Summary and Acknowledgments 

The CLEO Collaboration continues to produce results in charm physics using CLEO II, II.V, 
III statistics. We have a broad charm-physics dedicated program and are looking forward to 
CLEO-c/CESR-c data. We gratefully acknowledge the effort of the CESR staff in providing us 
with excellent luminosity and running conditions. This work was supported by the National 
Science Foundation, the U.S. Department of Energy, the Research Corporation, and the Texas 
Advanced Research Program. 
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